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Evolution of Galaxies and the Tully— Fisher Relation 

E. P. Kurbatov/fl A. V. Tutukov/fl and B. M. Shustov 1 '! 

1 Institute of Astronomy, Moscow, Russia 

We study the evolution of the [0/Fe]-[Fe/H] relation and the dependence of the 
iron abundance on distance from the galactic plane z in a one-zone model for a 
disk galaxy, starting from the beginning of star formation 1]. We obtain good 
agreement with the observational data, including, for the first time, agreement for 
the [Fe/H]-,z relation out to heights of 16 kpc. We also study the influence of the 
presence of dark matter in the galaxies on the star-formation rate. Comparison of the 
observed luminosity of the Galaxy with the model prediction places constraints on 
the fractional mass of dark matter, which cannot be much larger than the fractional 
mass of visible matter, at least within the assumed radius of the Galaxy, ~ 20 kpc. 
We studied the evolution of disk galaxies with various masses, which should obey 
the Tully-Fisher relation, M oc R 2 . The Tully-Fisher relation can be explained as 
a combination of a selection effect related to the observed surface brightnesses of 
galaxies with large radii and the conditions for the formation for elliptical galaxies. 



I. INTRODUCTION 

Many thousands of theoretical papers 
have been written on the evolution of galax- 
ies, often applying sophisticated mathemat- 
ical methods: iV-body models, multidi- 
mensional gas-dynamical models, statistical 
methods etc. However, a number of very im- 
portant results can be obtained using meth- 
ods that are relatively simple and are not 
computationally demanding. In [H 0, H El HI ; 
we suggested and implemented an approach 
that is both simple mathematically and rea- 
sonably self-consistent in the treatment of all 
complex physical processes determining the 
evolution of a galaxy. 

In mathematical terms, this evolution re- 
duces to the solution of two equations. One 
defines the star-formation rate (SFR) assum- 
ing the complete ionization of the gas com- 
ponent of the galaxy, which is taken to be 
homogeneously distributed over the galactic 
disk. The other equation describes variations 
of the thickness of the galactic disk based on 
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the condition of virial equilibrium, the in- 
put of mechanical energy into the gas com- 
ponent by super novae, and the dissipation of 
energy due to collisions of gas clouds. Nu- 
merical modeling must also take into account 
the evolution of the galactic stellar compo- 
nent: the return of matter to the interstellar 
medium by old stars after the formation of 
the endproducts of their evolution — black 
holes, neutron stars, and degenerate dwarfs 
- and the exchange of matter between the 
galaxy and intergalactic medium. 

Despite the simplicity of this approach, it 
led to a fairly complete picture of the evo- 
lution of disk galaxies — in particular, the 
Milky Way — that was consistent with ob- 
servations. This made it possible for the first 
time to distinguish the role of the loss of 
heavy elements from the Galaxy during the 
formation of the radial chemical-composition 
gradient, and to derive self-consistent distri- 
butions of metals in the z direction over a 
scale of ~ 2 kpc. This provided an explana- 
tion for the origin of enhanced metallicities 
(by a factor of a few over the solar value) of 
galaxies harboring quasars in their nuclei, via 
the enhanced density and higher _minimum 
mass of the stars in these regions 
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FIG. 1: [0/Fe]-[Fe/H] dependence for models with various initial accretion times r a . The solid 
curve shows the standard model (r a = 0). The filled circles show data on subgiants Q, open 
triangles data on subgiants and main-sequence stars 0], filled triangles data on sub dwarfs [j 1, 
inverted open triangles data on halo giants HI , and rectangles the positions of R CrB stars h], 11 ] . 
The remaining objects are HE 0107-5240 0], CS 22949-037 Q, CS 29498-043 Q, CS 29497-030 
Q, LP 625-44 Q, and G 77-61 Q. 



The present paper continues our study of 
the evolution of disk galaxies using the model 
developed by us earlier. This model enabled 
us to trace the chemical evolution of galaxies 
starting from their formation. However, the 
absence of data on very old, low-metallicity 
stars prevented comparison of the model re- 
sults with observational data for the initial 
stages of the evolution of the Galaxy. New 
data for stars with extremely low metallici- 
ties has recently become available, making it 
possible to trace the early stages of the en- 
richment of the Galaxy in heavy elements. 
Section 2 discusses the evolution of abun- 
dances implied by our theoretical modeling. 
Section 3 discusses modern data on the dis- 
tribution of metals in the direction of the 



Galactic rotational axis on scales of 16 kpc. 
In our previous studies, we used only fairly 
old data on stars that are located not further 
than 3 kpc from the disk — again hinder- 
ing application of the model to the earliest 
stages of Galactic evolution. Another impor- 
tant circumstance influencing the evolution 
of the Galaxy is the presence of dark mat- 
ter in the Galaxy. The nature and properties 
of dark matter are far from fully understood 
(see, e.g., the review [3]), but the necessity 
of taking this factor into account is beyond 
doubt. Section 4 presents an estimate of the 
fractional mass of dark matter obtained in 
our model. 

The good agreement of the model re- 
sults with the observations encourages us to 



3 



use the model to interpret various statistical 
relations (mass-metallicity, mass-luminosity, 
etc.). In recent years, possibilities for ap- 
plying the Tully-Fisher relation [T3|, which 
was discovered in 1977 by the two named 
authors, in various studies of galaxies have 
been widely discussed. This law relates the 
luminosity and rotational velocity of a disk 
galaxy. There are various ideas about the 
origin of the Tully-Fisher relation. In Sec- 
tion 5, we suggest an explanation that seems 
to us quite natural, based on an analysis of 
galactic evolution using our model. Finally, 
Section 6 presents a summary of our results. 

We made some modifications to the set of 
standard parameters for the Galaxy for the 
computations; these are described in the cor- 
responding Sections. 



II. EVOLUTION OF THE OXYGEN 
AND IRON ABUNDANCES 

Analyses of the oxygen and iron abun- 
dances are the classical methods for testing 
our understanding of the evolution of stars 
and galaxies. These elements were adopted 
as tools for directly analyzing the abundances 
of the products of type I supernovae (SN la, 
which produce mainly iron) and type II su- 
pernovae (SN II or SN Ib,c, which produce 
mainly oxygen). The time between the for- 
mation of a massive star and its explosion as 
a SN II or SN Ib,c does not exceed ~ 10 7 yr. 
According to one possible scenario, SN la re- 
sult from the mergers of degenerate dwarfs 
[lij ]. The lifetimes of SN la precursors cover 
a wide range, from 10 7 to 10 10 yr, with the 
most probable values being from 2 x 10 8 yr 
to 2 x 10 9 yr [19]. The possibility that SN la 
make a significant contribution to the pro- 
duction of iron in the Galaxy at ages lower 
than ~ 10 9 yr is supported by the discovery 
of an [O/Fe] ratio equal to the solar value 
in a distant quasar 20] with z = 6.4 (age 
~ 7 x 10 8 yr). We have used here a standard 
model for the evolution of the Galaxy: the 
closed model with outer radius R = 20 kpc 



and mass M = 2 x 1O 11 M described in jl|. 
The only difference is in the lifetime of SN la 
precursors (which was assumed to be 10 9 yr, 
close to the median value in the scenario for 
the evolution of close binaries [HI). In previ- 
ous computations of the standard model (H, 
this time was taken to be 3 x 10 8 yr. Test 
runs have shown that this change does not 
result in any fundamental changes in the re- 
sults of the computations; therefore, we shall 
adopt this modified model as our standard. 

We considered the evolution of the 
[0/Fe]-[Fe/H]_ relation in a number of pre- 
vious studies p, 0, E| , using stars in the solar 
neighborhood with relative iron abundances 
[Fe/H] > —3. The metallicity range consid- 
ered here is wider, extending to [Fe/H] = —7. 
The main results are presented in Fig. 
When [Fe/H] < —1, the abundances of oxy- 
gen and iron are determined primarily by 
SN II explosions, leading to an overabun- 
dance of oxygen by a factor of two to four. 
At an age of ~ 5 x 10 8 yr, SN la begin to ex- 
plode, producing primarily iron [2jJ (in the 
closed model, this age corresponds to [Fe/H] 
~ —0.5). The ratio [O/Fe] begins to de- 
cline with time (along with the increase of 
[Fe/H]), and approaches the solar value at 
an age of several Gyr. The model provides 
an acceptable distribution for most stars with 
[Fe/H] > -3. 

Among stars with the solar abundance 
of heavy elements, it has long been known 
that R CrB stars have abnormally high abun- 
dances of oxygen and carbon. The positions 
of R CrB stars with solar iron abundances 
are plotted in Fig. [TJ based on the data of 
Q; low-metallicity ([Fe/H] « -0.5- - 2.0) 
R CrB stars are plotted based on the data of 
[ill ] . The overabundance of oxygen in R CrB 
stars reaches nearly two orders of magnitude, 
and may be nearly independent of [Fe/H] 
(Fig. EJ. The reason for the overabundance 
of O and C in R CrB stars is the episodic 
penetration of the convective envelope of the 
star on the asymptotic giant branch onto the 
layer enriched in the products of helium burn- 
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ing |22J]. This process enhances the observed 
abundances of C, 0, and N in red super- 
giants by nearly two orders of magnitude. 
This oxygen-enrichment mechanism is possi- 
ble only for red (super)giants with luminosi- 
ties that are more than two orders of magni- 
tude higher than the solar luminosity. Stars 
of lower mass may become enriched in carbon 
or oxygen due to mass exchange between the 
components of close binaries. This may ex- 
plain the origin of such oxygen-rich stars as 
LP 625-44 and CS 29497-030 0. Stars with 
carbon abundances [C/Fe] = 2 and [Fe/H] 
= —4 are known l2i 



An increasing number of observations of 
stars with [Fe/H] < —3 have been reported 
in recent years. The very low metallicities 
of these stars reflect the earliest stages of 
the evolution of the Galaxy 24] . Such stars 
can be arbitrarily separated into two groups 
(Fig. [TJ. The first group contains stars that 
show a substantial excess of oxygen relative 
to the model: G 77-61 0, CS 29498-043 
0, CS 22949-037 0, CS 22957-027 
([O/Fel « 2, [Fe/H] = -3.11), and CS 31062- 
012 ([O/Fe] « 2, [Fe/H] = -2.55), and 
R CrB. The second group contains the giant 
HE 0107-5240 0. While the O excesses in 
stars of the first group can be explained by 
the penetration of their convective envelopes 
into the layer enriched in oxygen or by mass 
exchange between binary components, of the 
oxygen overabundance in the stars with the 
lowest metallicities may be primordial, in- 
dicating that they formed during the first 
~ 10 7 yr of the Galaxy's existence. Figure 1 
shows that the [O/Fe] ratio in the earliest 
stages of the Galaxy's evolution is very high, 
and decreases with age. Since SN la have not 
yet started to enrich the interstellar medium 
in iron, a natural question arises: What is 
the reason for this behaviour of [O/Fe]? 

The reason for the increase of the [O/Fe] 
ratio is the increase in the production of oxy- 
gen and the decrease in the production of iron 
with increasing initial masses of SN II pre- 
cursors. According to the numerical models 



|26l l27j , the amount of oxygen M Q produced 
by stars with masses of 15-120 M & can be 
written 



Mo 

M m 



0.01 



1.7 



where Mi is the initial mass of the star. The 
amount of iron produced by 13-25 M Q stars 
can be estimated as fej 



Mi 



Fc 



21 



Mi 



-1.87 



M© VMo, 
Combining these relations, we obtain 

Mo K 1n -tfMi 



Mi 



5 x 10" 



Fc 



Mr: 







3.57 



This expression is valid for 13 < Mi/M < 
120, since onl y s tars of these masses were 
studied in 0, 27]. Of course, when working 
with we must bear in mind that, while 
the production of oxygen by massive stars is 
estimated relatively accurately, the iron yield 
remains uncertain. It may depend on several 
parameters that are not accurately known, 
such as the initial chemical composition of 
the stars, the rotational velocity of the pre- 
supernova core, the magnetic-field strength 
etc. It is clear from observations that the 
iron yield may be even smaller than the mod- 
els |27J predict. For instance, SN 1997D, 
with a total mass of ~ 25M Q , ejected only 
2 x 10~ 3 M^ of nickel (which produces iron as 
it decays) [28j, much less than the ~ 0.05M Q 
expected from the models |27|. On the other 
hand, SN 2002ap, with an initial mass of 
~ 25M , produced 0.07M Q of nickel 0, 
close to the model expectations. 

According to (JTJ), obtaining the ratio 
[O/Fe] ~ 2 observed in the extreme low- 
metallicity star HE 0107-5240 (Fig. HJ) re- 
quires a supernova with an initial mass of 
~ 5OM and an age of 3 x 10 6 yr. Thus, judg- 
ing from the iron abundance, HE 0107-5240 
is among the first stars to be formed in the 
Galaxy, during the first several million years 
of its existence. Stars with [Fe/H] < —3 have 
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masses of the order of a solar mass. This 
suggests that the initial enrichment of the 
Galaxy in heavy elements was accomplished 
by ordinary stars with masses from ~ 1M 
to - 1OOM . 

The metallicity of the first Galactic stars 
could grow due to the accretion of interstellar 
gas. If the radius for capturing interstellar 
gas by a star is r = 2GMv~ 2 , where G is 
the gravitational constant, M the mass of the 
star, and v its relative velocity, the star can 
accrete during the lifetime of the Galaxy a 
mass 



AM 

Mo 



10" 



; n H 



30 



where nn is the current number density of the 
hydrogen in the gaseous phase (in cm -3 ) and 
v is the relative velocity (in units of 30 km/s). 
According to (J2J), even in stars that are com- 
pletely devoid of heavy elements, the rela- 
tive iron abundance could become equal to 
that observed in HE 0107-5240 (Fig. EJ) due 
to accretion (if there is a negligible role of 
the stellar wind, low mixing efficiency in the 
outer layers of the star, etc.). Note that 
low-metallicity stars could have initially be- 
longed to a former low-mass satellite that 
merged with the Galaxy in the past. The 
evolution of the chemical composition of low- 
mass spheroidal galaxies may end in an early 
phase, after loss of the gaseous component. 
Naturally, stars in such a galaxy could retain 
their high [O/Fe] values. 

For completeness, we note an additional 
possible means of enrichment of population- 
Ill low-mass stars (which are initially devoid 
of heavy elements) in metals. Most stars are 
formed in stellar clusters, which, as a rule, 
are disrupted almost immediately after their 
formation because of the loss of the gaseous 
component due to the formation of HII re- 
gions and supernova explosions |3l( . The first 
stars may also have formed in clusters. Then, 
some fraction of the lost gas, along with the 
products of the first, most massive, super- 
novae, would be captured by low- mass stars 



with the initial low-metal chemical composi- 
tion. The capture radius is r « 2GM*v~ 2 , 
where M* w M is the mass of the cluster 
star and v ex the velocity of the expanding 
gaseous envelope of the cluster. The fraction 
of the captured matter will then be 



a 



G 2 Ml 
vlR 2 



N' 



(3) 



where R is the cluster radius, % the free-fall 
velocity at the edge of the cluster, and iV 
the number of stars in the cluster (assuming 
that all have the solar mass). If we assume 
that v s = v ex , N = 10 3 and M Fe 0.01 M 
(M « 50 M 0)), then a solar-mass 
star can capture ~ 1O _8 M of iron. This 
produces a star with an iron abundance of 
[Fe/H] pa —5, similar to that observed in 
HE 0107-5240 (Fig. Q). It is obvious that 
a high oxygen abundance will be observed, 
since stars with masses exceeding ~ 5OM 
have [O/Fe] ~ 2, as was shown above. Note 
that not all the parameters in (JHJ) are well 
enough known to enable confident conclu- 
sions about the role of this mechanism in the 
enrichment of population-Ill stars in heavy 
elements. The accretion of interstellar gas 
enriched in heavy elements produced by the 
first stars may substantially complicate, or 
even make impossible, the identification of 
low-mass population-Ill stars. 

Note that the [O/Fe]- [Fe/H] evolutionary 
curve obtained for the closed model shown 
in Fig. H] represents an upper bound for the 
oxygen abundance. In low-mass disk galax- 
ies, the products of SN II explosions can 
efficiently leave the parent galaxies via the 
formation of supershells j32|, thereby reduc- 
ing the growth rate of the oxygen abundance 
in these galaxies. Since type la supernovae 
do not produce supershells, the iron abun- 
dance in low-mass disk galaxies increases at 
the usual pace. The example provided by 
the galaxy IZw 18 j^, which has active star 
formation with [Fe/H] = —1.76 ± 0.12 and 
[O/Fe] = -0.3 ± 0.3, shows that the loss 
of the products of type II supernovae may, 
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FIG. 2: Distribution of iron abundance with height above the Galactic plane in the standard model 
(solid curve) and in a model taking into account the loss of heavy elements due to SN II explosions, 
for fractional masses of heavy elements lost equal to 0.16 (dashed curve) and 0.5 (dotted curve). 
The solid circles show the data from 301. 



indeed, be very important for dwarf galax- 
ies. As our model computations show, al- 
lowing for the loss of type I supernova prod- 
ucts in massive galaxies does not lead to any 
substantial differences in the [0/Fe]-[Fe/H] 
distribution. Moreover, special runs of the 
model taking into account the finite time for 
the accumulation of the Galactic matter show 
that increasing the Galaxy's formation time 
substantially reduces the [O/Fe] ratio at low 
[Fe/H] (Fig. HJ. Increasing the amount of 
data on low-metallicity stars in the future 
may enable estimation of the time for the ac- 
cumulation of the Galaxy's mass. 



III. DISTRIBUTION OF HEAVY 
ELEMENTS ABOVE THE GALACTIC 
DISK 

Let us now consider the distribution of 
heavy elements with height z above the 
Galactic disk. Various observations show 
that the metallicity gradient with z mea- 
sured using field stars may reach — 0.65/kpc 
(according to data for the thin disk |34j ) ) . 
while the lowest value obtained from obser- 
vations of open clusters is — 0.34/kpc |35^ . 
At the same time, some observations of field 
and open-cluster stars show no evidence for a 
metallicity gradient 0, |3?J. To summarize, 
at z < 4 kpc, the gradient is in the range 
(—0.55 . . . — 0.33)/kpc, while it is in the range 
(-0.18. . .0)/kpc at larger z Q. 

In the standard model, we adopted the 
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TABLE I: Evolution of the main parameters of the standard Galaxy model with age (SFR is the star 
formation rate in M & /yr, H g the semi-thickness of the gaseous disk, M g the mass of the gas, Z the 
fractional mass of heavy elements, Td us t the optical depth of the dust, L opt = £boi Adust (1 — e _Tdust ) 
the optical luminosity with allowance for absorption by dust) 



t, 


yr 


SFR 


Hg, 


pc 


M g /M 


z 


[0/H] 


[Fe/H] 


Tdust 


-^opt 


/L & 


Lbal/L® 


3.6 


• 10 6 


34.3 


2.0 • 


10 4 


2 





10 11 


7.0 


• 10" 


-7 


-4.4 


-6.9 


0.0012 


1.1 • 


10 10 


1.1 


10 10 


4.8 


• 10 6 


34.4 


2.0 • 


10 4 


2 





10 11 


2.0 


• 10" 


-6 


-3.9 


-5.8 


0.0034 


1.2 • 


10 10 


1.2 


10 10 


6.3 


• 10 6 


34.5 
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10 4 
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4.3 


• io- 
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-3.6 


-5.1 


0.0072 
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10 4 
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7.5 
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10 4 
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-3.2 


-4.3 


0.020 
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10 10 


1.4 


10 10 


1.5 
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35.3 
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10 4 
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-5 
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-3.8 
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1.6 


10 10 


2.5 
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36.2 


1.9 • 


10 4 
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-2.7 


-3.3 
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1.7 


10 10 


3.4 


•10 7 


37.0 
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10 4 
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-2.5 


-3.1 
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• 10 7 


37.9 
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10 4 
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-3.0 


0.14 
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10 10 
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10 4 
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10 4 
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-4 


-2.2 


-2.7 


0.25 


2.0 • 


10 10 


2.3 


10 10 


9.8 


•10 7 


43.8 


1.5 • 


10 4 
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10 11 


2.2 


• io- 


-4 


-2.0 


-2.5 


0.36 


2.2 • 


10 10 


2.6 


10 10 


1.2 


• 10 8 


47.0 


1.4- 


10 4 
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10 11 


3.0 


• io- 


-4 


-1.9 


-2.4 


0.47 


2.3 • 


10 10 


2.9 


10 10 


1.7 


■10 s 


54.6 


1.2 • 


10 4 


1 


9 


10 11 


4.4 


• io- 


-4 


-1.7 


-2.2 


0.72 


2.5 • 


10 10 


3.5 


10 10 


2.5 


• 10 s 


69.7 


8.7- 


10 3 


1 


9 


10 11 


7.2 


• io- 


-4 


-1.5 


-2.0 


1.1 


2.7- 


10 10 


4.5 


10 10 


3.0 


• 10 s 


83.3 
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10 3 


1 


9 


10 11 


9.5 


• io- 
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-1.8 


1.5 
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10 10 


4.0 


•10 8 


123.7 
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10 3 


1 


8 


10 11 


1.6 


• io- 


-3 


-1.2 


-1.6 


2.3 


3.2 • 


10 10 


8.1 


10 10 


5.6 


■10 8 


243.5 


1.6 • 


10 3 


1 


5 


10 11 


3.1 


• io- 


-3 


-0.9 


-1.3 


4.0 


3.6 • 


10 10 


1.5 


10 11 


7.1 


•10 8 


342.7 


6.2 • 


10 2 


1 


1 


10 11 


6.9 


• io- 


-3 


-0.5 


-1.0 


6.5 


3.8 • 


10 10 


2.5 


10 11 


9.2 


•10 8 


161.7 


5.5 • 


10 2 


7.2 


10 10 


1.4 


• io- 


-2 


-0.2 


-0.6 


8.6 


2.4 • 


10 10 


2.0 


10 11 


1.2 


•10 9 


76.1 


4.9 • 


10 2 


4.7 


10 10 


2.1 


• io- 


-2 


0.0 


-0.4 


8.2 


1.6 • 


10 10 


1.3 


10 11 


1.6 


■10 9 


33.7 


5.4- 


10 2 


3.3 


10 10 


2.7 


• io- 


-2 


0.0 


0.0 


7.4 


1.2 • 


10 10 


8.9 


10 10 


2.1 


■10 9 


21.1 


5.0 • 


10 2 


2.5 


10 10 


3.3 
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-2 


0.1 
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10 9 


6.7 


10 10 


2.9 


■10 9 


12.8 


4.1 • 


10 2 


1 


8 


10 10 


3.8 
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-2 


0.1 


0.4 


5.6 


9.5 


10 9 


5.3 


10 10 


3.8 


•10 9 


7.2 


3.8 • 


10 2 


1.3 


10 10 


4.2 


• 10" 


-2 


0.2 


0.4 


4.4 


9.4 


10 9 


4.2 


10 10 



value of [Fe/H] clS db function of the disk 
semi-thickness for the dependence of the iron 
abundance on the height above the Galactic 
plane. This dependence is shown in Fig. El 
along with data from [3(j • We can distinguish 
three regions with different metallicity gradi- 
ents: d[Fe/B]/dz = — 3/kpc for z < 0.7 kpc, 
d[Fe/R]/dz = -0.2/kpc for 0.7 < z < 5 kpc, 
and d[Fe/R)/dz = -0.1/kpc for 5 < z < 
16 kpc. This corresponds to different iron 
abundances in stars of the thin disk, thick 



disk, and halo, and also to different spatial 
scales within these components. 



Computations of galactic evolution allow- 
ing for the loss of heavy elements (Fig. |2J) and 
accretion have shown that the "openness" of 
the galaxy weakly influences the model de- 
pendence of [Fe/H] on z. Thus, our one-zone 
model can explain well the variations of the 
distribution of heavy elements with z. 



8 



IV. EVOLUTION OF A DISK 
GALAXY: THE INFLUENCE OF 
DARK MATTER 

The main results obtained for the stan- 
dard model are listed in the table. Com- 
parison with observations indicates that this 
model provides an adequate description of 
the chemical and dynamical evolution of the 
Galaxy, as well as of its star-formation his- 
tory. For example, the data of 38] pro- 
vide evidence that young galaxies with M* ~ 
10-100M o /yr have optical depths of ~ 5-10 
in the direction perpendicular to the plane 
of the disk, as is confirmed by our model. 
The large optical depths of distant, and thus 
young, galaxies could influence the extra- 
galactic distance scale (which is based on 
the assumption that SN la are standard can- 
dles), and therefore conclusions concerning 
the regime of the cosmological expansion. We 
should also bear this in mind in connection 
with analyses of the brightness of supernovae 
at z = 2-5. Variation of the time to accu- 
mulate the mass of the Galaxy in the initial 
stage of its formation did not influence sub- 
stantially the Galaxy's parameters after the 
Hubble time (Fig.QJ. 

Our model has another important param- 
eter — the mass of dark matter, i.e., of mat- 
ter (baryonic or otherwise) that cannot be 
directly observed, but provides a (sometimes 
dominant) contribution to the gravitational 
field. The fractional mass of dark matter 
is small on small scales. For instance, esti- 
mates of the fractional mass of dark matter 
in the Sun yield less than 2-5% 39] , and this 
is probably true for other stars as well. On 
scales of 100 kpc or more, the fractional mass 
of dark matter may be much larger [lfil ]. 

We studied the role of dark matter in the 
evolution of the Galaxy by changing the ex- 
pression for the gravitational potential in the 
standard model to 

M dm 



E 



GM g H 2 (M c 



R 2 



+ 



R 



the mass of gas, H g the semi-thickness of the 
gaseous disk, Mq the mass of visible mat- 
ter, the semi-thickness of the stellar disk, 
Md m the mass of dark matter, and R the ra- 
dius of the disk. Figure El shows the evolu- 
tion of the main parameters of the Galaxy as 
functions of the ratio of the masses of dark 
and visible matter, which was varied from 
zero (corresponding to the standard model) 
to two. Reduction of the mass of matter in- 
volved in star formation results in a delay 
of the star-formation burst, and so reduces 
the thickness of the gaseous disk. This is 
in qualitative agreement with estimates fol- 
lowing from the observational relation be- 
tween the relative thickness of the disk and 
the mass-luminosity relation for galaxies |4£j . 
Increasing the mass-to-luminosity ratio in the 
model decreases the galaxy's relative thick- 
ness. The observed optical luminosity of the 
Galaxy places constraints on the fractional 
mass of dark matter, < 50%. Note that 
our initial mass function, dN oc M~ 2,35 dM 
for M > 0.1 M & , overestimates the num- 
ber of low-mass stars compared to the ob- 
served value. According to j4l|, the observed 
mass function changes slope at M ~ 1M Q 
(dN oc M~ h5 dM for M < 1M Q ), which 
makes the fraction of stars with < 1M less 
pronounced. Thus, an excess of matter ap- 
pears in our model compared to the obser- 
vations. A simple analytical estimate of its 
fractional mass is ~ 37%. 



V. TULLY FISHER RELATION 



About a quarter of a century ago, Tully 
and Fisher discovered a dependence be- 
tween the luminosity of a disk galaxy and its 
rotational velocity. In modern form, this re- 
lation can be written 



AM B = (8±l)Alog« 



rot j 



(4) 



where G is the gravitational constant, M g 



where Mb is the absolute B magnitude of 
the galaxy and v TOt is the rotational velocity 
of the galactic disk, estimated from the flat 
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FIG. 3: Evolution of the Galaxy with allowance for dark matter. SFR is the star-formation rate, 
Mg/Mc the ratio of the masses of gas and visible matter, H g the semi-thickness of the gaseous 
disk, Z the heavy-element abundance, Td us t the optical depth of the dust, and L opt the optical 
luminosity with allowance for absorption by dust. Computational results are shown for ratios of 
the masses of dark and visible matter M^/Mq equal to (solid curve), 0.5 (long dashed curve), 
1.0 (short-dashed curve), and 2.0 (dotted curve). 
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FIG. 4: Mass-radius diagram for various galaxies; galaxies obey the Tully-Fisher relation in the re- 
gion T-F. The lines tsf = Tsn and tsf = THubbie define the positions of galaxies with star-formation 
timescales tsf and th, respectively. The region LSBG is occupied by low-surface-brightness galaxies 
i3 . Early-type galaxies are located in the region EG 44, 4fJ (see Sect. 5). 



part of the rotation curve. We can use the 
mass-luminosity relation for disk galaxies, 
M/M Q ~ 1OL/L 0, to estimate the rela- 
tion between the mass of a disk galaxy M and 
its radius R. Assuming v rot oc M 1//2 i? _1 / 2 , we 
obtain 

R oc M°- 5±a07 . (5) 

Relations similar to (jlj) for varioius groups of 
disk galaxies have been studied by numerous 
authors (see, e. g., 0, EH 4t|). If selection 
effects are not considered, it appears that the 
locations and slopes of the log R-logM re- 
lation derived from (@J) are fairly uncertain 
(they vary by about a factor of three for 
galaxies with masses from 10 7 to 10 n M Q ). 
The Tully-Fisher relation can be written 



M = ttTjR 2 , where E is the surface density 
of gravitating matter. The observed range 
of S is rather wide. According to [42], disk 
galaxies display values S ~ 10-600 M Q /pc 2 ; 
another estimate yields £ = 2O-15OOM /pc 2 
|50| . Figure 0] plots a mass-radius diagram 
for galaxies of various types that obey the 
Tully-Fisher relation (located inside the "T- 
F" region in this figure), elliptical and early- 
type galaxies (EG; |44l |45| ) . and low-surface- 
density galaxies (LSBG~0, 0). 

We used (0) to compute a series of mod- 
els with masses of 10 6 - 4.5 x 1O U M . The 
results show that, in this range, the mass of 
a galaxy exerts almost no influence on the 
star-formation history or the evolution of the 
chemical composition. In particular, the fi- 
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nal heavy-element abundance does not de- 
pend on the mass, in contradiction with the 
observed increase of the heavy-element abun- 
dance with increasing galactic mass and lu- 
minosity This may be explained by the 
loss of heavy elements by low-mass galaxies. 

Our formalism 0, Hf enables us to derive 
a simple relation defining the time scale for 
star formation in a spherical galaxy 



r SF = 3 x 10 s 



R\ 



M- 



ii 



where R\ is the galactic radius in units of 
10 4 pc and M\\ the mass of its gaseous com- 
ponent in units of 10 11 M Q . The observed 
masses and radii of young galaxies are con- 
sistent with r SF = 10 7 -10 9 yr The 
spherical model describes the early stages of 
the evolution of disk galaxies. We can place 
limits on the time for the existence of star- 
forming galaxies. An upper bound for the 
region they can occupy in Fig. 4 is given by 
the Hubble time: r SF = r H = 1.4 x 10 10 yr, 
as is shown in Fig. 0J In galaxies located 
above this line, the SFR is so low that only 
a small fraction of the gas has been turned 
into stars over the Hubble time. Note that 
the position of this bound virtually coincides 
with the position occupied by low-surface- 
brightness galaxies ji^, El • Nearby galaxies 
with masses of 10 6 -10 12 M Q from the catalog 
|H3 | abut on the upper boundary of the region 
of star-forming galaxies (they are not shown 
in Fig. E}; t ne star- formation time scale in 
such galaxies is rs F ~ 80 9 -10 10 yr. 

On the other hand, the SFR in high- 
density galaxies is so high that almost all 
the gas has been transformed into stars even 
before the first explosions of SN II, whose 
precursors have lifetimes tsn ~ 5 x 10 6 yr. 
This corresponds to the relation r SN < r SF , 
or Mn > 60i?4 (the line r SN = r SF is also 
plotted in Fig. HJ). As soon as the SN II begin 
to explode in these galaxies, they clean out 
any remaining gas. This is the scenario for 
the evolution for elliptical galaxies. The mass 
lost by old stars in these galaxies is unable to 



revive star formation, since SN la explosions 
drive the galactic wind, hindering star forma- 
tion throughout the galaxy, with the possible 
exception of a circumnuclear region with a 
high gas density (see, for example, [EH). In 
fact, only a few percent of SO and E galaxies 
show signs of star formation in their nuclei 
[56(. Some SO and E galaxies may form via 
collisions of disk galaxies 57, 58[. However, 
the number of such galaxies is not large, since 
observations 59 show only a small increase 
in the fraction of E galaxies with the age of 
Universe; hence, most formed at large red- 
shifts. Albeit with a fairly large scatter, the 
positions of early-type galaxies in Fig. HI are 
close to the boundary defined by the rela- 
tion r SF = r SN . However, elliptical galaxies 
formed in collisions of disk galaxies may also 
have larger sizes. 

Figure 0] shows that the galaxies with ac- 
tive continuing star formation are located in a 
strip defined by the conditions tsn = rg F and 
tsn = th; the slope of this strip is close to 
the slope of the Tully-Fisher relation. This 
suggests that the power of the R oc M a law 
is determined by the similarity of the slopes 
of the boundaries of the region occupied by 
disk galaxies in the mass-radius diagram. At 
the lower boundary of this zone, galaxies that 
are initially elliptical are formed, while the 
upper boundary is apparently defined by two 
circumstances. The first is the condition that 
there be a sufficient SFR on a time scale that 
is shorter than the Hubble time. The second 
is the near coincidence of the upper boundary 
with the position of galaxies that have surface 
brightnesses at the detection limit. Galaxies 
with lower surface densities, S < 10M Q /pc 2 , 
are probably not detectable. 



VI. CONCLUSION 

Our study of the evolution of a disk 
galaxy with star formation governed by ion- 
ization has enabled us to apply this model 
to several new fields. In particular, we have 
shown that high oxygen-to-iron abundance 
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ratios are characteristic of the very first stars 
formed: [O/Fe] « 2 when [Fe/H] < -5. This 
is a result of the reduction in the iron produc- 
tion and increase in the oxygen production 
with increasing initial mass of type II super- 
novae. This suggests that some stars with low 
iron abundances and relatively high oxygen 
abundances (G 77-61, HE 0107-5240, Fig. CJ) 
may be either the earliest second-generation 
stars or first-generation stars that are "con- 
taminated" by products of the first super- 
novae. Our model distribution of metals with 
height z above the Galactic disk agrees with 
the observed distribution to z ~ 16 kpc. In 
our model, the well-known Tully-Fisher rela- 
tion may stem from a combination of obser- 



vational selection effects and conditions that 
are necessary for the formation of disk galax- 
ies. 
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